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Abstract 

An analysis of the role of general relativistic effects on the decay of neutron 
star's magnetic held is presented. At hrst, a generalized induction equation 
on an arbitrary static background geometry has been derived and, secondly, 
by a combination of analytical and numerical techniques, a comparison of 
the time scales for the decay of an initial dipole magnetic held in flat and 
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curved spacetime is discussed. For the case of very simple neutron star mod- 
els, rotation not accounted for and in the absence of cooling effects, we find 
that the inclusion of general relativistic effects result, on the average, in an 
enlargement of the decay time of the field in comparison to the flat space- 
time case. Via numerical techniques we show that, the enlargement factor 
depends upon the dimensionless compactness ratio e = , and for e in the 
range (0.3 , 0.5), corresponding to compactness ratio of realistic neutron star 
models, this factor is between 1.2 to 1.3. The present analysis shows that 
general relativistic effects on the magnetic field decay ought to be examined 
more carefully than hitherto. A brief discussion of our findings on the impact 
of neutron stars physics is also presented. 
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I. INTRODUCTION 



It is well known |IJ that a magnetic field in a plasma of finite conductivity is subject to 



diffusion and dissipation. Diffusion leads to a spreading of inhomogeneities while dissipation 
is due to the Ohmic decay of the currents producing the field. More concretely, a magnetic 



Accordingly, if L is a typical length scale of the field structure, then it will decay-diffuse 
in a characteristic time scale rohm given by: rohm = ^ c t ■ Depending upon the prevailing 
conditions, the Ohmic decay time r 0hm can range from seconds, in the case of a copper 
sphere of radius of a few centimeters up to rohm = 10 10 years or even much longer for 
astrophysical settings, as in the case of the sun |j] or a neutron star 0. 

The interactions of large scale cosmic magnetic fields with plasmas is a problem of great 
importance in astrophysics and cosmology. A particularly thorny issue nowadays concerns 
the origin and maintenance of cosmic magnetic fields. Although large scale fields have been 
observed |3|], a satisfactory explanation of their origin is still lacking. Peebles considers the 
issue of the origin of the primordial magnetic field as one of the most important unsolved 
problems in cosmology. At the same time the gigantic field of the pulsars begs for an 
explanation || . The general consensus of the astrophysical community [fj is that such large 
scale fields have been generated via an episode of dynamo action [0, and then gradually 
suffer Ohmic decay due to the finite conductivity of the medium. It appears therefore that 
an understanding of the factors influencing the decay of large scale fields, combined with 
relevant observations, may offers important clues towards a better understanding of the 
initial scale involved as well as clues regarding its origin. 

In neutron stars the decay of the magnetic field is an issue of out most importance by 
itself J8| and accordingly there has been an intense effort by astrophysicists to understand 



field B(t,x) in a plasma of uniform conductivity a evolves, in flat space-time, according to 



the following diffusion equation M: 
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the factors governing this decay. As far as we are aware all theoretical modeling of magnetic 
field decay in neutron stars utilized the familiar flat space time form of Maxwell's equations 
(an exception to this rule constitutes the recent work of ref. ||). Although the employment 
of such framework is a fruitful one and provides us with valuable informations, it altogether 
neglects the background curvature of the spacetime which for the case of neutron stars is 
not any longer weak. It would be worth to stress in that regard that curvature can modify 
considerably flat space-time solutions of Maxwell's equations. For instance, the reader may 
compare the solution describing a dipole magnetic field on a Schwar zschild background 
|T0| to that of a flat space time. The presence of the logarithmic term in the former (see 



eqs. |3.17| further below) is a sole consequence of the non vanishing curvature. This example 
suggests that the role of the spacetime curvature on the decay process of magnetic fields 
ought to be examined more thoroughly than hitherto. In that respect, we are aware only of 
the recent work of Sengupta, [[|, where an investigation of general relativistic effects in the 
magnetic field decay of neutron stars have been attempted. However this work is restricted 
to the study of magnetic fields confined only to the outermost layers of a neutron star and 
furthermore it is assumed that those outermost layers (and thus also the magnetic field B), 
are embedded on a Schwar zschild background geometry. Thus strictly the framework of ref. 
H, deals exclusively with magnetic decay on a Schwar zschild background. In addition to 
those approximations and according to the sentence following Eq. 15 of Sengupta's second 
work, the author fails to include general relativistic effects on the outer boundary condition 
for matching the inner field with the outer vacuum dipolar field across the surface of the star. 
In contrast, in the present work, a broad framework dealing with general relativistic effects 
on the magnetic field decay on an arbitrary static geometry and with proper allowance of 
the correct general relativistic inner and outer boundary conditions is presented. Moreover, 
and in contrast to the approach of ref. |J, we formulate the entire problem avoiding the 
introduction of a vector potential and the associated ambiguities. Our analysis shows that 



general relativistic effects (TJJ can influence the field decay, but the precise manner that 
this influence manifests itself depends upon the class of observers called in to describe the 



field decay. For the magnetic field of a non rotating neutron star it is natural to describe 
the field decay relative to the class observers that find themselves at rest relative to the 
star, ie the class of Killing observers. Relative to such observers, we find that relativistic 
effects are influencing the field decay via two major modes: the gravitational red shift as 
well as the intrinsic curved geometry of the spatial sections constituting the rest space of 
the Killing observers. Subsequent numerical analysis shows that the red shift factor is the 
dominant one in slowing down the field decay. Overall we find that the inclusion of relativistic 
effects make the decay time of the field larger than, but of the same order of magnitude, 
as in flat space-time. Nevertheless the preliminary study of the present paper utilizing a 
simple non rotating neutron star models suggests that general relativistic effects should 
be given further considerations. We explicitly illustrate the impact of relativistic effects 
upon the magnetic field decay, by examining the evolution of a magnetic field permeating a 
constant density neutron star, first in their presence and secondly without them. Although 
for both treatments we have obtained exponential decays, the decay time in the presence of 
relativistic effects, on the average, is enlarged by a factor that depends crucially upon the 
value of the compactness ratio e = ^rjf- Specifically for values of e in the domain (0.3 , 0.5), 
characterizing realistic neutron star models, we find that the decay time is 1.2 to 1.3 larger 
than the corresponding flat decay time, while for higher values of e, it can be larger. We 
may add parenthetically that the term "average" increase in the decay time, is explained in 
details in section (IV) of the paper. 

The present paper is organized as follows: In the following section, starting from 
Maxwell's equations on a static spacetime we first derive the relevant induction equation 
taking into account the curved nature of the background spacetime geometry. It should be 
stressed however that the employment of a static geometry does not leave room for incor- 
porating gravitomagnetic (Lense-Thirring) effects in the induction equation, as the latter 
would manifest themselves relative to non static backgrounds, but we do hope to present 
such analysis in a future work. In section III, we specialize the induction equation to a 
simple neutron star model and a detailed analysis of the content of the induction equation is 
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presented. In the same section the sensitive issue of the boundary conditions accompanying 
the induction equation is also addressed. In the section IV, we discuss numerical solutions 
of the curved spacetime induction equation and an assessment of the relativistic factors 
influencing the field decay is discussed. Furthermore in the same section, a comparison of 
the field decay in curved and flat spacetime is also presented. In the concluding section, a 
brief discussion of the physical implications of our results to neutron stars physics is pre- 
sented and possible extension of the present work is outlined. Finally we have included an 
Appendix where a few intermediate calculations leading to the main equations of section II 
are presented. 



II. INDUCTION EQUATION ON A STATIC BACKGROUND GEOMETRY 



Maxwell's equations, in covariant form, are as follows [I2",|T3|: 

47T 



V a F r 



a/3 



J, 



ft 



(2.1a) 







(2.1b) 



where F a p = —Fp a , J a and V are the coordinate components of the Maxwell tensor, the 
conserved four current and the derivative operator respectively. Given a solution F a p of the 
above eqs., an observer with four velocity U a , U a U a = — 1, measures electric and magnetic 
fields (E, B) with corresponding coordinate components given respectively by: 



F a pUP 



^ s F l5 U? 



(2.2) 



where e Q/ g 7 5 stands for the four-dimensional Levi-Civita tensor density fI3 |. We shall be con- 



cerned in the present paper with particular solutions of |2.1| where the current J is described 
by the following relativistic extension of Ohm's law, as it was first formulated by Weyl ||15|| : 

J a = ag^F^V 1 



(2.3) 



where in above equation, (V a , a) stand for the four velocity of a conducting neutral plasma 
and its scalar electrical conductivity |16| respectively. Although eqs. [O] to |2.3| are valid for 
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any kind of background geometries and plasmas characterized by arbitrary four velocity and 
conductivity, hereafter we shall restrict our consideration to background geometries that 
are globally static. Staticity in turn allows us to select coordinates so that the spacetime 



geometry can be written in the form (see for instance discussion in ref. [p~2] , p~3[1 ) : 



ds 2 = -e 2q, (dx°) 2 + ^dx l dx 3 (2.4) 

where x° = ct, 7^ are functions of the the spatial coordinates x l , (i = 1,2,3), and by £ 
denote the hypersurface orthogonal timelike Killing vector field obeying: £ a £ a = — e 2 *. For 



the above form of the line element, Maxwell's equations |2.1| and the current conservation law 
V a J a = can be re-written in an equivalent form involving only the components (E l , B l ) 
of the electric and magnetic fields respectively, as well as the charge density cp = —U^J^ 
and spatial current density J 1 as measured by the Killing observers fT7Ul8|| . More precisely 
if by we denote their four velocity then eqs. |2.1| yield the following equivalent set (see 
Appendix for details, or ref. |I7||L8||): 

DiE? = 4vrp , = (2.5a) 

e^D 3 (ZB k ) = -Zr + — (2.5b) 

e^ Dj (ZE k ) = -— (2.5c) 

jjh^P)_ + D .ji + ji D .i ogZ = (2.5d) 
ox^ 

where in above D stands for the covariant derivative operator associated with 7, e 1 ^ represent 
the (coordinate) components of the three dimensional totally antisymmetric Levi-Civita 
tensor density defined on the x° = const slices and Z = (— £ a £ a )^ = e* is the red shift factor 
which in the language of in the 3 + 1 approach to spacetime or (and) electrodynamics, is 
also refered as the lapse function fTBfl . 

With Maxwell's eqs. in the above form we can derive an induction equation by repeating 
the same steps leading to the derivation of its flat counterpart (see for example discussion 



in J!]]). For a plasma at rest relative to the Killing observers, combined with Ohm's law 



and the MHD approximation (i.e. neglecting the displacement current |0| from the right 
hand side of |2.5b| ), one obtains from |2.5[ -abc the following form of the generalized induction 
equation: 



dx° J 



Aira 



e k lm D x {ZB r , 



(2.6) 



This last equation describes the time evolution of a magnetic field configuration that 
finds itself in a conducting medium. In principle one could write down the explicit form 
of the dynamical evolution equation once a choice of background geometry has been made. 
However before we do so, we would like to make a further specialization of the eqs. [2.5| and 
[2.6| , so that their interrelationship to the familiar flat space three plus one formalism of 
Maxwell's eqs. is more transparent. Here, following the spirit of |18[ and particularly ||20| . 



we shall sacrifice the manifest three-covariance of eqs. [2.5| and [2.6| with respect to arbitrary 
coordinate transformations of the t = const sections, for the benefits of practical usefulness. 
As was pointed out in ref. [ |TB| , |2T)| ] , if one defines suitably the components of (E,B) and 



under some weak constraints upon the spacetime geometry, then Maxwell's equations can 
be recast in a more "user friendly" form. This new form employs concepts familiar from 
the language of the three dimensional vector analysis expressed in orthogonal curvilinear 
coordinates and such approach to curved spacetime electrodynamics is particularly useful 
for astrophysical purposes. Having in mind further astrophysical applications of our results 



we shall recast eqs. |2.5| -abc in such a form. Such form requires that the geometry of the 
spacetime permits the introduction of coordinates so that the spatial three element ds 2 3 ^ of 
[2.4] could be recast in the following form: 

ds 2 {3) = h\(dx 1 ) 2 + h\{dx 2 ) 2 + h 2 3 (dx 3 ) 2 (2.7) 

where the scale factors hi = hi(x , x 2 , x 3 ) are for the moment arbitrary functions of 
(x 1 ,^ 2 ,^ 3 ). In the Appendix, (see also |TE| , |2T)|| ) we show that for such geometries eqs. |2T3 



can be written in the following form: 
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V-£ = 47rp, V-S = 



(2.8a) 



Vx(ZB) 



4tt z j IdE 
c c dt 



(2.8b) 



(2.8c) 



V- J + J-V{logZ) = 



(2.8d) 



where we have written the current conservation law for an electrically neutral plasma and 
in above equations the symbols (V- , Vx , V) stand for the divergence, curl and gradient 
operators respectively, expressed entirely in terms of the scale factors hi (see Appendix for 
their explicit representation). We also remind the reader that all vector components in 
equations |2l] are physical frame components taken with respect to the field of orthonormal 
frames = ^-^i, (i — 1,2,3), naturally singled out by the line element |2.7| . Using now 
eqs. |2.8|, or directly from eq. |2.6|, upon eliminating the coordinate components of B in favor 



of its frame components, the induction equation takes the following form: 

ldB 



c dt 



Vx 



Ana 



Vx(ZB) 







(2.9) 



Equations |2.8| and |2.9| are the main equations of this section. In the special case of a 
Schwar zschild background, naturally they are reduced to those of ref. [ 2"0[] , and in the case 



of a plasma of uniform conductivity the generalized induction eq. |2.9| reduces to eq. [0| in 
the limit of flat space. 



III. MAGNETIC FIELD DECAY INTERIOR TO NEUTRON STARS 

In the neutron star's interiors the MHD approximation is well justified [^] and we shall 
explore the content of the relativistic induction equation |2.9| , by applying it to study the 
evolution of magnetic fields associated with neutron stars. Since the main purpose of the 
present work is to investigate the impact of the spacetime curvature upon the magnetic field 
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decay, as a first preliminary step we shall adopt a rather simplified neutron star model. The 
chosen model primarily avoids technicalities that may obscure the issue at hand but at the 
same time it shows clearly the potential impact of the curvature on the magnetic field decay. 
Accordingly, and to avoid laborious numerical computations, we shall ignore the rotation 
of the neutron star and thus shall adopt as the background geometry a non-singular, static 



and spherically symmetric one. Hence, the scale factors of eq. |2.7| will be taken as: 



= / _ 2G^y = ( 1 _?my 1 , hl = s , K = r > ^ s (3 ,, 



while for the moment the lapse or red shift factor Z = Z(r) = e*^ and the "mass function" 
m = m(r) are arbitrary functions of the radial coordinate. 

We shall begin our analysis of the magnetic field decay by assuming that at some initial 
time t a an axially symmetric distribution of a magnetic field B(t Q , r, 9) permeates the entire 
star. We are not concerned here upon the mechanism that brought such a field into existence 
but rather we are interested in its evolution. Its evolution is considerably affected by the 
electrical conductivity a, but as a part of the adopted simplified picture and in order to 
emphasize the effects of space-time curvature we shall take a to be spherically symmetric 
and shall ignore any cooling effects that may influence its temporal evolution. For an axially 
symmetric field B, it is convenient to decompose it into the so called poloidal Bm and 
toroidal part Buy In terms of the orthonormal basis vectors (e r , eg, e^) those parts are 
defined respectively by: Bu^ = B r e r + B eg and B^ = B^e^ with (B r ,B 6 ,B^) arbitrary 
functions of (t, r, 9) respectively. One can then easily conclude from the induction eq. |2.9| 
that, as long as the scalar conductivity is spherically symmetric, the toroidal and poloidal 
parts of B, evolve independently of each other f21f . Such decoupling is rather convenient 
since it implies that if the initial distribution of the magnetic field is purely poloidal then 
it will not develop a toroidal component in the course of its evolution and vice versa. For 
simplicity, in the present paper we shall examine the effects of the spacetime curvature only 
on the evolution of a purely poloidal field Bu^ = B r e r +B e eg. For such field B, it follows from 
2.8b| that the current J is along the direction, and thus the current conservation eq. [2.8d 
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is identically satisfied. Furthermore via Ohm's law, and use of |2.8b| (with the displacement 
current ignored), it follows that the electric field E is a purely toroidal and axisymmetric 
field and Gauss law V-E = is satisfied as well. Consequently, from the system of eqs. |2.8| , 
we are left to satisfy the constraint V--B = 0, solutions of which will be evolved by the 
induction eq. 

Taking into account the poloidal and axisymmetric nature of B as well as the formula 



of the div operator V-, listed in the Appendix, in view of the scale factors of [3.1|, one easily 

finds that V-B = implies: 

(l _ 2M Y 1 9(r 2 B r ) + 1 d(B e sm9) = Q 
\ r J r dr sin# 89 

We shall look for separable solutions of the above equations in the form: 

B r = F{t, r)Q x {6), B e = G(t, r)6 2 (0) (3.3) 

with the functions F, G, 6i, 9 2 to be determined. Substituting the above representations 
of (B r , B 6 ) in |3T^ and separating variables one gets the following equivalent system: 

2MV1 9 (^)_ AG = 



r dr 
1 <9(sin#0 2 



+ A9i = (3.4b) 



sin 6 89 

where A stands for a separation constant. The second equation can be solved in terms of 
the Legendre polynomials by taking A = 1(1 + 1), I = 0, 1, 2..., and: 

2 = sinfl^P^ Qi = -Pi(y), y = cos9 (3.5a) 
dy 



On the other hand, for such A, eq. |3.4a| is satisfied provided, for 1^0, one chooses G(r,t) in 
the following form; 

1 / 2M\ 5 1 d(r 2 F) 

G ' i,r ' = ((iTT) ;V < 3 ' 5b > 

We shall disregard the / = mode since, as it is clear form above, it corresponds to a 
monopole field B. With the exclusion of monopole fields, the components of an arbitrary 
axisymmetric poloidal field can be written as a superposition of " 1-poles" in the form: 
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(3.6a) 



V ' ' ; ^«=i Z(Z + 1) V r / r <9r rfy 



(3.6b) 



To simplify algebra, and on physical grounds, we shall restrict our considerations to the 
detailed analysis of only the I = 1 mode. Such mode corresponds to a dipole field and 
such configuration is expected to be present and dominant within neutron stars. For 1=1, 
eqs. |3]6| yield: 



1 / 2M\ ^ d(r 2 F) 

B r (t, r, 9) = -Fit, r) cos 9, B 6 \t,r,9) = — 1 ^— ^ sin 9 

2r V r J or 



(3.7) 



where for notational simplicity we write here after F instead of F\. On the other hand, for 



any poloidal axisymmetric field with components (B r ,B ), the induction equation |2.9| on 
the background geometry of [2.7|, yields the following two non-trivial evolution equations: 



dB r 1 
+ 



d 



dx° hohs 89 



cA 
Ana 







(3.8a) 



dB e 



1 d 



dx° h r h$ dr 



cA 



(3.8b) 



where: 



A 



h„ 



h r hg 



d_ 

Or 



(h e B e z) 



d_ 

89 



{h r B r Z) 



(3.8c) 



When one now inserts in eq. |3.8a| the explicit forms of the components of (B r , B ) corre- 



sponding to a dipole field in the form shown in eq. |3.7| , as well as the scale factors of eq. | 3~T 
then gets: 



4vr(T dF 

c dx° 



1 - 



2M\2 1 d 



-) 



r 2 dr 



Z 1 



2M\5 d(r 2 F) 



dr 



2ZF 



(3.9) 



In arriving at the above equation we have taken explicitly into account the spherically 
symmetric nature of the scalar conductivity o. We may point out that for non spherical a, 



the right hand side of 3.9 contains gradients of a along the meridian directions but for our 



12 



simple neutron star model a spherical conductivity is rather adequate. On the other hand 
identical manipulations of (3.8b) leads to: 

2 dF c 



d_ 

Or 



dx° Ana 



2My d 
r J dr 



Z 1 



2My d{r 2 F) 
r ) dr 



+ 2 



Ana 



-ZF 



(3.10a) 



from which we infer that 



dF 



dX° 4:71(7 



2M\2 d 
Or 



Z 1 



2M\ 2 d{r 2 F) 
Or 



+ 2 ZF = g(9, (ft, t) (3.10b) 



4.7TCT 



where g(8,<j),t) an integration "constant". A comparison then between (3.9) and (3.10b) 
shows that necessary g = 0. If we further define: F = r 2 F, then one finds either from |3.9| 



or 



3. 10b| that F satisfies: 

47T(T dF 



c dx° 



2My d 
r J dr 



Z 1 



2M\2 dF 
dr 



(3.11) 



The above equation essentially describes the evolution of the the dipole field components 



22] as they are measured by the Killing observers. Linearity of Maxwell's and induction 



equation implies that |3.11| specifies a unique solution up to an arbitrary rescaling. This 
rescaling freedom will be fixed later on by a suitable matching of the interior dipole field to 
a corresponding asymptotically vanishing exterior dipole one. 

Taking Z — 1, M = in eq. |3.11| one recovers the standard equation describing the 
evolution of a dipole axisymmetric poloidal field in flat space namely: 

AnadS d 2 S 2S 



c 2 dt dr 2 



(3.12) 



where in order to avoid confusion we have denoted the analogue of F(t, r) for the flat space 
time case by S(t, r). The latter function often in the astrophysics literature is referred to as 
the Stoke's function |24[]. In order to get some insights into the significance of the various 
terms appearing in (3.11), the general relativistic counterpart of |3.12| , we shall rewrite the 
former in an equivalent form so that a clean comparison between the two could be afforded. 
Eliminating the (areal) radial coordinate r in favor of the physical proper radius l(r) of the 
r = constant spheres, via dl = dr(l — — )~a J Eq. [3.11| takes then the following form: 
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4:71(7 dF d OF 2Z ~ 

= ~ iw ( a) 

A comparison then to eq. |3.12| shows that relative to the Killing observers, general rela- 
tivistic effects can influence B-decay in three ways. Namely via the presence of the red shift 
factor Z, its gradient and as well as via the intrinsically curved nature of the rest space of 



the Killing observers, ie the t = const hyperfaces. The latter manifest itself in |3.13a| via 
the term r(l), a term which in general satisfies r{l)^l, implying that that the rest spaces 
of the Killing observers are intrinsically curved. From the above mentioned three factors 
the spatial gradient of Z makes negligible contribution to the field decay and this has been 
verified numerically. Neglecting this gradient then equation |3.13a| takes the following form: 

Ana dF _ d 2 F 2F 

In this form a clean comparison to equation (3.12) can be afforded. The right hand sides 
of the two equations involve physical spatial gradients and the differ only by terms of the 
order 0(^^). On the other hand their left hand sides as they stand cannot be compared. If 
however one reasonably replaces Z(l) by some averaged value < Z >, then the left hand side 
of (3.13) involves also physical temporal gradients. In that event one gets a first flavor of the 
magnitude of the general relativistic effects. They modify the corresponding flat spacetime 
results by terms of order unity. Of course such conclusions has to be also documented at 
the solution level as well, and as we shall further ahead this indeed is the case. 

Having thus identified the manner by which relativistic gravity effects the magnetic field 
decay, our assignment is now to access the relative importance of each of the above two 
factors. In the following section we shall do so by resorting to numerical computations. 
However, before we pass to that issue let us first record the suitable boundary conditions 
to be imposed upon the corresponding S(t, r) , F(t, r) in order to describe sensible physics. 
The required conditions for both, ie the flat and the general relativistic case, are drawn 
by demanding that the interior B ought to be a non singular field at all times and at all 
spatial points, and in addition it ought to join smoothly across the surface of the star to an 
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exterior asymptotically vanishing dipole field. For the flat space case, taking M = and 
Fit, r) = — |f in eq. |3.7| one gets: B = cos 9 e r — sin 9 ee from which we infer that the 
magnitude B 2 of the interior magnetic field is given by: \B\ 2 = ^ cos 2 9 + ^-(|^) 2 sin 2 9. 
Accordingly a regular field at the star's center requires lim s ^ to be finite as the center 
of the star is approached. On the other hand an asymptotically vanishing dipole magnetic 
field in flat space due to a magnetic moment \i, is described by [p]]: B = 2 ^° s 6 e r + MS1 ° eg. 
It follows then from the above expressions that a C° matching of the interior magnetic field 
to slow varying exterior dipole field f25|, requires that across the star's surface, ie at radius 



R, S(t,r) should satisfy: 

= -S (3.14) 



R — 

dr 



R 

The relatively simple nature of (3.12) as well the simple form of the boundary- regularity 
conditions outlined above, permit us to construct exact closed form solutions. In fact, it is 
not difficult to verify that for the case of a star of a uniform conductivity cr, a sequence of 
exact solutions of eq. (3.12) obeying the above described conditions is given by PB |: 



S n it, x) 



sin n« cos mre 



n 2 7r 2 x nn 



e T " = ipn(x) e t " (3.15) 



where r n = Au ? 2 = \ ZQ ifS x — % and n takes the values (1, 2, 3....). 

The above sequence of exact solutions offers a clear picture regarding the behavior of a 
magnetic field in a conducting medium that finds itself in a flat space time. Constructing for 
instance the field Bi(t, r, 9) corresponding to Si(t, r), one immediately sees that an observer 
at fixed (r, 9) finds that the magnitude of Bi(t, r, 9) decays exponentially with a characteristic 
e-folding time given by T\ = ^§-- On the other hand expanding an arbitrary initial field 
configuration B(to,r,9) in terms of the eigenf unctions [ip n ,n = 1,2..), one can easily see 
the spatial diffusion of the initial distribution. For a plasma characterized by an arbitrary 
a, although the decay and diffusive nature of the initial B field remain intact, it is rather 
difficult to estimate analytically the characteristic decay time as well as to find out whether 
the decay will be channeled into an exponential phase. It is sufficient, however, to stress 
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that as long as we are in flat spacetime the decay process is controlled by the conducting 
properties of the background medium and, of course, the length scale of the initial field 
distribution. 

Let us now turn the discussion to the formulation of the appropriate conditions to be 
imposed on solutions of eq. |3.11| . Since as already indicated in the introduction section, a 
dipole field on a Schwar zschild is modified considerably from its flat form and, as eq. |3.11| 
shows, relativistic effects modify the local behavior of the relativistic Stoke's function F(t, r), 
one expects modification of the boundary conditions as well. As far as the behavior of F(t, r) 
at the star's center is concerned, by arguing in the same manner as in the flat space case, a 
non singular dipole field requires F(t, r) to satisfy identical conditions at the star's center as 
its flat counterpart, namely fi'm ^v should be finite as the center of the star is approached 
(after all, the principle of equivalence holds). Although this is the case at the star's center 
the corresponding boundary conditions at the star's surface are markedly different. Recalling 
that the frame component of the vector potential A = A^ = A^e^ describing a magnetic 
dipole on a Schwar zschild background is described by [ |10| : 



A„ 



3/i sin^ 

AM 2 



x In 1 



1 



x 



1 

2x 



+ 1 



(3.16) 



where x 



2M(R) 



R <r < oo, and fi is the dipole moment. From |3.16j , one then obtains the 



corresponding field B = B r e r + B e cq where the physical components (B r , B ) as measured 
by the Killing observers are given by: 



2fi cos 6 



3x 3 m(l - x' 1 ) + 3a; 2 + -x 



(3.17a) 



B h 



/ism 9 



6x 3 (l - oT 1 )^ ln(l - aT 1 ) + 6x 2 



1-i 



x 



(3.17b) 



A comparison of (3.17, ab)) with the corresponding eq. |3.7| and a C° matching between the 
two along the star's surface, requires that the exterior dipole magnetic moment /x should be 
identified with the generally slow time varying part of the function F(t,r) Moreover 
the gradient F(t, r) along the radial direction should obey: 
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dF(t,r) 



dr 



= G(y)F(t,R) (3.18a) 



with: 



g(y) = y I , y~ "A (s.isb) 

|/ 2 ln(l -y !) + y + 2 



aIld ^ 2MWY 

Thus the behavior of the interior dipole field in the presence of curvature is described by 
F(t, r) satisfying the differential eq. p.ll| , subject to boundedness of F ^ as the star's center 
is approached, and additionally obeying |3.18a at its surface. Before we turn our discussion to 



the construction of solutions of eq. |3.11| subject to the above discussed conditions, we would 
like to write an explicit formula for the time evolution of F(t, r) under the assumption that 
geometry of the spacetime corresponds to a static spherically symmetric star, solution of 
Einstein's equations. We may recall that in the derivation of eqs. (3.11) and (3.18,ab) we 
assumed an arbitrary non-singular, static, spherically symmetric background geometry with 
the only constraint that it joins smoothly to an exterior Schwarzs child field across the 
surface of the star. Nowhere in the derivation we needed the explicit form of the M(r) = 
nor the form of Z = Z(f). Hereafter we shall become more explicit and shall take the 
background interior geometry to be a non singular solution of the coupled Einstein-perfect 
fluid system. As is well known, and under the assumption that the background Maxwell field 
makes negligible contribution to the structure of the star [27|, static spherically symmetric, 



perfect fluid solutions of Einstein's equations imply satisfaction of the following differential 
equations between the metric functions $(r), M(r), the hydrostatic pressure P(r), and mass 
density p{r) (see for instance fl2|]l3|l ): 



rf$ _ M(r) + 4irr 3 P(r) 

~dr~ ~ r 2 (i _ m^i j " (3 - 19a) 

4nr 2 p (3.19b) 



dr 

dP(r) . . Mir) + Anr^P 

^ = -^ +j V(i_m^ < 3 - 19c) 
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Making use of those equations, and restoring the fundamental units, we obtain from |3.11 
the following equation to be satisfied by the relativistic Stokes function: 



4™ e _ Hr) dF_ = f 1 _ 2Gm{r) \ d 2 F | 1 
c 2 dt \ c 2 r I dr 2 r 2 



2Gm(r) 4ttG , ( P(r) 



cr \ cr 



(3.20) 



where for typographical convenience we shall write here after F(r,t) instead of F(r,t). The 
above equation via |3"77| , describes the evolution of any axisymmetric, dipole, poloidal field B 
that finds itself interior to a spherical perfect fluid star. In the above form it includes all three 
relativistic factors influencing the field decay. Since the distributions of (m(r), P(r), p(r)) 
are related via Einstein's equations directly to the spacetime curvature, eq. |3.20| shows 
implicitly that the influence of spacetime curvature on the decay of the magnetic field is a real 
effect and cannot be removed via coordinate transformations. In principle one could insert 
in p.20| the appropriate distributions of m(r), P(r), p(r) resulting from integrating the 
Oppenheimer-Tolman-Volkov equation, specify o = a(r,t), and construct the history of the 
5-decay. We shall report elsewhere our findings of this rather laborious numerical integration 



For the purpose of the present paper we shall integrate |3.20| for a rather simple system 
introduced and discussed in the following section. The goal of this section is to show that 
the general relativistic eq. |3.20| (or its approximate forms corresponding to eq. (3.13b) 



under the assumption of a uniform conductivity admits decay modes analogous to those of 
the flat space case with one important difference: The corresponding e-fold decaying times 
are longer in the relativistic case. We interpret this amplification of the e-fold decay times 
as resulting from the non vanishing space time curvature. Unfortunately however it is not 



easy to construct analytically the exact decaying modes of the full relativistic system |3.20 
or its approximate versions (3.13a), and thus we shall resort to numerical computations. 
The emphasis in those computations is the probing of the dependence of the corresponding 
e-folding times upon the value of the red shift factor or (and) upon the strength of the 
curvature of the spatial sections. 
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IV. MAGNETIC FIELD DECAY IN A CONSTANT DENSITY STAR. EXPLICIT 

RESULTS. 



We shall consider in this section the decay of a magnetic field in a neutron star of constant 
density. The assumption of a constant density star, although not a very reliable approxima- 
tion of a real neutron star, offers the advantage that the Einsteins equations can be solved 



analytically (see for instance [[T^,|T3|]), and thus provides us with closed form expressions for 
the coefficients of the induction equation, eq. |3.20| , and the boundary condition, eq. |3.18| . 
In particular, in Box 23.2 of the ref. |13| the distribution of the various hydrodynamical 
and geometrical variables are plotted as functions of the areal coordinate r. As we have 
already discussed, the exact decaying modes of the flat space-time induction equation for a 
uniform conductivity are explicitly known (given by eq. |3.15|) , while the corresponding de- 
caying modes of the full curved space-time equation |3.20| are presently unknown. We shall 



therefore resort to numerical techniques in an attempt to get insights into behavior of the 
space of solutions of eq. |3.20| . 

Viewed as an initial-boundary value problem, (3.20) is a diffusive initial value problem for 
which the standard numerical technique is the Crank-Nicholson implicit integration scheme 
(see for example ref. [29] for a description). We checked our numerical code by evolving the 



fundamental mode of the flat space-time case (ie take n=l in eq. |3.15| ) and compare the 
numerical solution with the analytical one: we obtained an accuracy better than 1% until 
times up to 10 • rj, with n = = Z£ ^f I the corresponding decay time of the n = 1 flat 
fundamental mode. 

Before we turn our discussion of the numerical results it is helpful to view the time 
evolution of a chosen initial distribution from a complementary point of view. On general 
grounds F(t = 0, r) as well as its time evolution can be (formally) expanded in a series of 
the following form: 

c 2 X n t 

F(t,x) =Y,a n e-^?g n {x) (4.1) 
where the summation is extended over all eigenmodes g n (x) of the corresponding (singu- 
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lar) Sturm-Liouville eigenvalue problem arising from eq. [3.20| and the associated boundary- 
regularity conditions: 



Lg n + \ n e g n = 



(4.2a) 



1 - 



2Gm(x)\ d 2 



c 2 x 



dx 2 



^ X 2 



2Gm{x) AnG 3 f P(x) N 

C 2 C 2 I C 2 ^ i 



dx x 2 (4-2b) 



where in the present case the coefficients in L are determined by the geometrical and hydro- 
dynamical variables of the constant density star solution, x = and R is the areal radius of 
the star. It follows now from (4.1) that if the eigenvalues are positive and well spaced, then 
after t » where Ai is the lowest eigenvalue of the above system, then the evolution 
of the distribution will channel into an exponentially decreasing phase with the dominant 
contribution in the sum (4.1) coming from the "first" term. Our subsequent described nu- 
merical computations exhibits such feature and this property allows to construct numerically 
the lowest eigenvalue of the above system |SU . 



In the following numerical calculations we have taken the areal radius to be R = 10 
km, a constant uniform conductivity o = 10 25 s _1 typical of neutron star values (which 
implies rohm = 4.44 10 9 yrs), and we consider various neutron star masses characterized 
by different values of the dimensionless compactness ratio: e = ^jj^f = 0, 0.3, 0.4, 0.5, 0.6, 
0.740, 0.810, 0.865, and 0.889. The first one corresponds to a flat background space time, 



the last four are those values used in the numerical plots of ref. |[L3|| , while current realistic 
neutron star models are characterized by e in the range 0.3 to 0.5 [^TJ. For each value of 
e, at first we have solved numerically the full relativistic induction equation |3.20| , by taking 
the initial F(t = 0,r) to be equal to the Stoke's function S(t = 0,r) of the corresponding 
first fundamental decay modes of |3.15| (ie take n = 1 and t = in |3.15| ) . After performing a 
long time-integration of eq. |3.20| subject to the conditions cited earlier on, we find that the 
evolution of F(t, x) channels into an exponentially decaying mode which means, according 



to 4.1, that the evolution of the initial distribution eventually is described by the first non 
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vanishing term in the series expansion |4.1| . This behavior of F(t, x) allows us to determine 
only the lowest eigenvalue Ai of (4.2) from our numerical outputs. Besides the explicit 
determination of Ai, our numerical treatment allows us to construct the magnetic field as 
well. In Fig. [1], we plot as a function of coordinate time t, the magnetic field as perceived 
by a Killing observer located at the star's pole for the various values of the compactness 
ratio. Fig.(l) shows that once curvature effects are incorporated and upon ignoring the 
initial transit time during which the field is in a superposition of various curved eigenmodes, 
the field follows an exponential decay law (as would have done in the absence of gravity) 
but now the corresponding e-folding time is longer than the corresponding flat spacetime 
case. Thus even though we have started with identically prepared systems their evolution 
is distinct, a distinction traced in the influence of relativistic effects. It should be stressed 
however that the content of Fig.l does not by itself provide us with a clear overall picture of 
the field decay. It rather provides us with a characteristic physical decay time as perceived 
by a Killing observer situated at the surface of the star and this decay time should not 



be extrapolated as being the physical decay time over the entire star [R2j. In fact, each 



Killing observer located at some r, will compute a physical e-fold decay time r(r) given by: 
r(r) = Z(r)(Ai) _1 = Z ^ r ^ hm and obviously this value changes across the star. (In this 
formula, we have parametrized Ai so that (3 = 1 corresponds to flat spacetime). Because of 
this spatial dependence of r(r), in order to get a better insights into the dynamics of the 

2 

decay, in Fig.2, we have plotted Ai = j3— — , as a function of the compactness ratio e. In 
the same figure, for comparison purposes, we have plotted the value of the red shift factor 
at the stars center [Z = e®^) and surface (Z s = e®^) respectively. Thus it follows from 
Fig. (2), that the relativistic corrections to the lowest eigenvalue Ai, are bounded from above 
by Z s while from bellow (almost) by Z . It is more instructive however and complements the 
content of Fig. (2), a plot showing the physical decay times as measured by Killing observers 
located at the center and at the surface of the star respectively. Fig. (3) stands for such plot, 
and its content shows that the physical decay time can vary considerably across the star. In 
fact there are regions around the star's center, where the physical decay time is shorter than 
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the corresponding flat space case and this effect is more pronounced as the compactness ratio 
increases. In contrast to what occurs in the vicinity of the star's center, in the crust region 
the physical decay is always larger than the corresponding flat case. Because of this behavior, 
largely due to gravitational time dilation effect, we assign an overall physical decay time, 
by averaging the physical e-fold decay at the center and the surface of the star respectively. 
This amounts to assigning an overral a red shift factor Z for the entire star equal roughly 
to its value at the middle of the star. With this type of averaging, Figs. (2,3) shows that for 
small values of the compactness ratio the overall physical decay time is almost identical to 
the flat space time case. However, as the compactness ratio increases, the relativistic effects 
become more apparent. For the case of neutrons stars with range in the realistic domain 
ie e in the range (0.3 , 0.5), and via the averaging procedure outlined above, the overall 
physical e-fold decay time is (1.2 — 1.3) larger than the corresponding flat case. Although 
the content of Fig. [L], Fig. |^ and Fig. || show the impact of relativistic effects upon the field 
decay, by themselves they do not offer a clear insight as which (if any) of the two factors, 
ie red shift or spatial curvature, are responsible for the dominant contribution in the field 
decay. In order to access their relative importance we solve numerically eq. |3.13a| in two 



extreme cases [^3J and show the numerical outputs in Fig. (4). First eq. |3.13a| is solved under 
the assumption r(l) = I and in this approximation the relativistic effects on the decay are 
solely due to the red shift factor Z(l). In fig. (4) the numerical outputs are indicated by the 



label: "curved time". In the opposite extreme, we adopt Z(l) = 1 in eq. |3.13a| and take r(l) 
as given by the metric corresponding to a constant density star. Thus in this approximation, 
the only relativistic effect influencing the decay is due to the spatial curvature. The resulting 
numerical outputs in Fig. (4) are marked by the label " curved space" . It follows then clearly 
from the content of Fig. (4) that for a constant density neutron star, the dominant effect in 
the field decay is due to the red-shift factor Z , since the corresponding eigenvalues indicated 
by "curved-time" graph are much closer to the corresponding exact eigenvalue indicated 
by "curved space-time" in Fig. (4). Moreover, the dominance of the red-shift factor holds 
through for all values of the compactness ratio e, and increases as e increases. 
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From the analysis presented so far it is clear that the more compact the star is, the longer 
is the e-folding time. As a consequence one expects that models of pulsars with soft equation 
of state to maintain a strong magnetic field for longer period of time than the corresponding 
models with a stiff equation of state. In turn such slow B field decay implies additional 
source of heating ie Joule heating and such additional heating, may explain the relatively 
high temperature observed in old neutron stars. However based on the present analysis, it 
is rather premature to draw definite conclusions. For instance cooling effects leading to the 
temporal variation of the conductivity as well as the the detailed structure of the star and 
its rotation has to be taken into account. Such study is currently under way and we expect 
to report in a future communication. 



V. CONCLUSIONS 

The behavior of the surface magnetic fields of neutron stars is a complicated and con- 
troversial issue. Many processes are believed to influence its magnitude and its subsequent 
evolution. Trapping for instance, of the field in the superconducting core is one possibility. 
The expulsion of the field out of this region is a delicate matter involving many different 
branches of physics |34]. Another possibility that in principle influences enormously the 



magnetic properties of neutron stars is related to the accretion processes immediately after 



the core collapse [35]. Accretion and particularly hypercritical accretion, can submerge the 
field of the new born neutron star beneath a layer of accreting matter thus in principle pro- 
ducing a delayed switched on mechanism for the pulsar activity ]36| . Furthermore according 



to recent work the neutron star may never turn on as a pulsar |37| if the accretion is hyper- 
ctical. Besides the above mechanisms influencing the evolution of neutron star's magnetic 
fields, many more have been introduced and discussed at length in the current literature. 
In this work, we have present a limited framework taking into account the effects of space 
time curvature on the field decay. For the simple neutron star models with a corresponding 
compactness ratio in the range (0.3 , 0.5), considered in the present work, we have seen an 
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overall increase in the decay time, (1.2 to 1.3) times larger than the flat spacetime value. 
Although the present work is preliminary and to assess the new effect more work is needed 
|28|j , it point towards to the direction that in a strongly gravitating system, effects due to 
space time curvature should not be neglected. 
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APPENDIX A: (3 + 1) FORM OF MAXWELL'S EQS. ON STATIC SPACETIMES 

In this Appendix we shall sketch a derivation of eqs.(2.8) starting from the covariant 



form of Maxwell's eqs. [2.1| . The derivation makes use of the existence of the hypersurface 
orthogonal timelike Killing field and although all the following computations can be done in 
a covariant fashion ]18|,^0j for brevity we work explicitly in the coordinate gauge of 2A . We 



shall also present formulas required for the derivation of equations in section II. 

Starting from the temporal component of the inhomogeneous Maxwell eq. [2.1a| , combined 
with the line element [2.4| and taking into account the fact that = —Fi°e® one immediately 
obtains: 

D i E i = - — J°e* = 4vrp (Al) 
c 

where we have defined the charge density p measured by a Killing observer by: cp = —J^U^ 
and we denote by D the covariant derivative operator associated with the Riemannian metric 
of the t = const spaces. Due to the fact that the Maxwell tensor F a p admits the following 
easy verifiable decomposition: F af3 = U a Ep — UpE a + e^sU 1 B s one gets the following 

24 



expression for its spatial part Fif Fy = e ijiB l U° = eijiB 1 . Passing now to the spatial 
components of p. la] one gets 



f) Air 

— -&D s {ZBi) = --ZJ* (A2) 

where we have introduced the red shift factor Z via Z = (— £ a £ a ) 5 — ^ instead of e*. On 
the other hand the second pair of Maxwell's equations |2.1b| can be written equivalently as: 



dx x ~^ dx^ ~^ dx v ^ ^ ^ 

Taking now all the indices to be spatial, and eliminating F^ one gets: DiB 1 = The other 
information encoded in the second pair of Maxwell eqs. can be revealed by considering the 
following arrangement of the spacetime indices: (/i, u,X — m, n, x°). For such arrangement 
one obtains: 

d(e mnl B l ) d(-U E n ) d(U Q E m ) = Q 

dx° dx m dx n ^ ' 

from which one easily obtains: 

3d! 

— + e lmn D m (ZE n ) = (A5) 
ox° 

The current conservation eq. V M J M = after a trivial rearrangement yields: 



To pass into the equivalent set (|2.8| -abcd) and eq. ^]9| involving physical orthonormal 



components, we project all tensors involved onto the natural set of orthonormal vectors (e^) 



and one forms (e l ), (i = 1,2,3) respectively, associated with the line element [2/?]. Thus 



for instance the electric field E can be written as: E = E 1 -^ = E l ei where E % = hiE 1 
expresses the relationship between coordinate and frame components of the vector field E 
and it is understood that no summation is involved over the repeated indices. With the help 
of the orthonormal components one for instance may rewrite eq. |A1| in terms of orthonormal 
components. Writing 7^ = /ii/i2^3 and eliminating the coordinate components in terms of 



the frame components of E, equation [Al| takes the following form: 
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where V- stands for the familiar divergence operator expressed in arbitrary orthogonal curvi- 
linear coordinates defined by the line element [2.7| . Similarly DiB 1 = can be written as 
V-B = 0. As far as the other set of Maxwell's eqs. are concerned, one can proceed in a sim- 
ilar manner. For instance starting from A2, one first multiplies the corresponding equation 
by the scale factor hi thus leading into: 

d{E l ) d{h k B~ k Z) An ■ 

-d^~h x h 2 h, dxi ~~~ JZ (A8) 

Recalling that the orthonormal components of the curl operator of an arbitrary three 
dimensional differentiable vector field A are given by ||38|| : 

one is lead immediately into eq. |2.8b| used in the text. Note also that the action of the 
gradient operator V acting on scalars is defined via: 

v/ = L2L + L2L + L2L (mo) 

hi dx 1 h 2 dx 2 h 3 dx 3 

Also in deriving eqs. (2.8-abc) of the main text we have used the following properties of the 
unit basis vectors (e*): e^-e^ = <5«, e^xe-- = e^e^ and the normalization e-^ = 1. We may 
also indicate that for typographical convenience the caret-symbol over frame components 
of the various tensors has been dropped. In particularly all vector, tensor components 
appearing anywhere in the main text after eq. |2.7| , are frame components. 
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[36] A. Muslimov and D. Page, Astrop. Journal, L77, 440, (1995), idid 458, 347, (1995) 
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FIGURES 



FIG. 1. Figure 1 

Dipolar field decay for a uniform density star in curved and flat space-time. The hor- 
izontal axis represents (coordinate) time in units of the flat space-time ohmic decay time 
T h m = 47iaR 2 /c 2 and the vertical axis shows the value of B/B = B(t,r = R,8 = 0)/B(t = 
0, r = R, 9 = 0). All models have the same areal radius (R = 10 km) and a constant uniform 
conductivity (cr = 10 25 s _1 ) typical of neutron star values (which implies T h m = 4.44 10 9 
yrs). The values of the compactness ratio 2GM/Rc 2 is indicated on each plot. The initial 
field profile is taken as the n = 1 eigenmode for the flat space-time, Eq. |3.15| , in all cases. 
The graphs show quite clearly the exponential decay in flat space-time, with r = T h m j-R 2 
within numerical accuracy, while, as expected, in curved space-time the decay initially devi- 
ates from an exponential law but rapidly converges toward the corresponding fundamental 
mode. 



FIG. 2. Figure 2 

The horizontal axis stands for the dimensionless comactenes ratio 2GM/Rc 2 , while the 
vertical axis corresponds to the values of (3. The graphs marked by e*° , e* s stand for the 
value of the red shift factor at the center and surface of the star respectively. 



FIG. 3. Figure 3 

The horizontal axis stands for the dimensionless compactenes ratio 2GM/Rc 2 , while the 
vertical axis corresponds to the ratio of the physical decay time r p h over the corresponding 
flat value r//. The graphs marked as "surface", "center" respectively, represents T ^? h — ^ 



Tfl P 

= ^jp while the corresponding horizontal line through (0,1) stands for the flat case 



FIG. 4. Figure 4 



34 



The horizontal axis stands for the dimensionless comactenes ratio 2GM/Rc 2 , while the 
vertical axis corresponds to the values of f3. The graphs marked as "curved space", "curved 
time" provides the eigenvalues corresponding the the case where Z = 1 and r(l) = I respec- 
tively, as explained in the text, while the graph marked as "curved space time" corresponds 
to the exact equation. 
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